Limits on primordial power spectrum resolution: An inflationary flow analysis 
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We investigate the ability of current CMB data to reliably constrain the form of the primordial 
power spectrum generated during inflation. We attempt to identify more exotic power spectra that 
yield equally good fits to the data as simple power-law spectra. In order to test a wide variety of 
spectral shapes, we combine the flow formalism, which is a method of stochastic model generation, 
with a numerical integration of the mode equations of quantum fluctuations. This allows us to 
handle inflation models that yield spectra that are not well described by the standard spectral 
parameterization. Using the latest WMAP data-set, we find a high degree of variation in possible 
spectral shapes. In particular, we find strongly running spectra arising from fast-rolling inflaton 
fields providing equally good fits to the data as power-law spectra arising from slowly-rolling fields. 
Current data poorly constrains the spectrum on scales k < 0.01 ZiMpc -1 , where the error due 
to cosmic variance is large. Among the statistically degenerate models, we identify spectra with 
strong running on these larger scales, but with reduced running at smaller scales. These models 
predict values for the tensor-to-scalar ratio, r, that lie outside the 2-a confidence interval obtained 
from SDSS+WMAP data for spectra that are parametrized as power-laws or spectra with constant 
running. By considering more generalized power spectra, we therefore open up regions of parameter 
space excluded for simpler models. 
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I. INTRODUCTION 



Understanding the nature of the cosmological den- 
sity perturbation is essential for the development of in- 
flationary cosmology. Recent years have seen a surge 
in the development of precision cosmology, with im- 
proved measurements of the temperature and polariza- 
tion anisotropics of the cosmic microwave background 
(CMB) as well as increasingly accurate large scale struc- 
ture (LSS) surveys. The task of faithfully reconstruct- 
ing the physics of inflation from this data is the cen- 
tral goal of inflationary phenomenology. The fundamen- 
tal physics responsible for inflation remains elusive, how- 
ever. Models of inflation abound with motivations stem- 
ming from supergravity, the string landscape, D-branes, 
extensions to the standard model and others p], 0, 
A reliable reconstruction of the inflaton potential, and 
therefore the physics responsible for inflation, hinges on 
accurate measurement and subsequent analysis of cos- 
mological data. There has been much progress in the 
development of methods of primordial pow er s pec trum 
reconstruction 0, i, 1, 0, S, I M EL M El, Hi SI , 
as well as methods for reconstructing the inflaton poten- 
tial [H, [H El, 0, [13Jpl- With the recent release of 
the WMAP3 data-set]22jand with future CMB missions 
eminent [H, [24[ , it is both important and timely to in- 
vestigate how well current and future experiments will 
be able to constrain the form of the primordial power 
spectrum. 

Recent analyses report a nearly scale invariant spec- 
trum of adiabatic, gaussian density perturbations, con- 
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sistent with the simp lest models of single field slow-roll 
inflation [H, [H, Hf [13, [H, HI • However, experimen- 
tal and theoretical inaccuracy, in addition to the limited 
detection of the polarization signal, limit parameter esti- 
mation efforts. The dominant source of error is that due 
to cosmic variance, a statistical uncertainty related to the 
fact that there is only one universe to observe. This effect 
is present in both temperature and polarization spectra 
and is most pronounced on large scales (at low CMB 
multipolcs). Cosmic variance results in a fundamental 
limit on the resolution of CMB spectra, which restricts 
the accuracy with which we can reconstruct the physics 
of the early universe. Since the primordial perturbation 
spectrum provides one of the few pieces of observational 
evidence for inflation, these uncertainties impose a limit 
on how well we can reconstruct the physics of the infla- 
tionary era. While single field slow-roll inflation is well 
supported by current data, it is of interest to determine 
whether more exotic inflationary scenarios might exist 
within the resolution limit of observational data. 

In this paper, we study our current ability to con- 
strain the form of the power spectrum. Using the 
latest WMAP data-set [12], we generate a best- fit 
base cosmological model characterized by 7 parameters, 
(ilih 2 , £l c h 2 , h, t, A, n s ,r). We are interested in identify- 
ing alternative models with more general power spectra 
that yield comparable fits to current data. Such statis- 
tically degenerate models can be considered effectively 
indistinguishable, and differences in the underlying infla- 
tion models unresolvablc. In order to test a wide range 
of spectral shapes, we introduce a Monte Carlo recon- 
struction prog ram that combines the inflationary flow 
technique 30, l3ll |32| with a numerical integration of 
the mode equations of quantum fluctuations. Aside from 
being a novel reconstruction technique in its own right, 
we find it highly suitable for the present investigation. 
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We use this method to stochastically generate an ensem- 
ble of power spectra to 6 t/l -order in the flow space. For 
each spectrum it is possible reconstruct the exact infla- 
ton potential. Because we solve for the spectrum nu- 
merically, we can successfully reconstruct spectra arising 
from slow-roll as well as strongly non-slow-roll inflation 
models. Additionally, the flow approach allows us to en- 
force that each model yield the appropriate amount of 
inflation, N g [46,60]. 

For each spectrum, we fix the non-spectral parameters 
at the best-fit values and calculate the new model's like- 
lihood. Wc establish a resolution criterion based on the 
relative likelihoods of the trial models and the best-fit 
model. By fixing the non-spectral parameters, we are 
not performing a Bayesian analysis of the enlarged pa- 
rameter space of flow parameters. By allowing all the 
parameters to vary as is done in Bayesian parameter es- 
timation, it is possible that new parameter degeneracies 
would allow for a wider range of possible spectral shapes 
to be in agreement with the data (see, for example, Ref. 
f3o|). Our method therefore only samples a subset of pos- 
sible spectra and should be viewed as a more conservative 
approach. 1 This study is more akin to a model selec- 
tion analysis: we attempt to identify models constructed 
from the inflationary flow space that are equally consis- 
tent with the data as models parameterized as simple 
power-laws. We utilize the p-value calculated from the 
model's chi-squared per degrees of freedom as a measure 
of statistical significance. This approach is frequentist in 
nature, yet we find general agreement with the Bayesian 
study of [13. 

This paper is organized as follows: in section II we re- 
view the inflationary flow formalism and in section III 
we introduce the new method which combines the flow 
approach with a numerical integration of the mode equa- 
tion. In section IV we investigate primordial power spec- 
trum resolution with the WMAP3 data-set. We find that 
with TT/TE/EE data alone, our ability to reconstruct 
the inflaton potential is strongly affected by the limited 
spectral resolution. We find that strongly running spec- 
tra arising from non-slow-roll inflation models are cur- 
rently indistinguishable from a best-fit power-law model 
typical of slow-roll inflation. However, since these non- 
slow-roll models are typically accompanied by large ten- 
sor components, any information regarding the B-mode 
polarization signal will significantly improve resolvability. 
Furthermore, we find that the established confidence lim- 
its on the tensor-to-scalar ratio obtained from SDSS and 
WMAP3 when running is included are violated by the 
non-power-law spectra identified in this study. We find 



spectra with values for r at k — 0.002 hMpc~ that lie 
outside the 2-cr contour of the n s -r marginalized likeli- 
hood [H, [2?| yield equally good fits to the data as the 
best-fit power-law. Section V includes discussion and 
conclusions. 



II. SINGLE FIELD INFLATION AND THE 
FLOW FORMALISM 

The evolution of a Fricdmann-Robcrtson- Walker 
(FRW) universe dominated by a single minimally- 
coupled scalar field (the inflaton) with potential V((f>) 
is given by the equations 



H = 



8tt 
3m P1 

8tt 
3m|j 



V{(t>) - ¥ 



(1) 

(2) 



We have assumed a flat universe described by the met- 
ric g^v = diag(l, — a 2 , — a 2 , — a 2 ), where a(t) is the scale 
factor and dots denote derivatives with respect to coor- 
dinate time, t. The evolution of the scalar field follows 
from stress-energy conservation and takes the form of a 
Klein-Gordon equation with a cosmological friction term, 



3i?0 + V' = 0, 



(3) 



where primes denote derivatives with respect to the field, 
4>. Equations ([T]) and can be combined to yield the 
alternative equations of motion, 
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where the Hubble parameter, written as a function of 
(j>, becomes the dynamical variable. The field <fi serves 
as a convenient time variable so long as it is monotonic. 
The second of the above two equations is known as the 
Hamilton- J acobi equation, and may be written more sim- 
ply as 



H 2 W 



1 



( 8tt 



\3mp[ 



where the parameter e is defined as 



'pi 



4tt \H(4>) 



(6) 



(7) 



1 A Bayesian parameter estimation analysis was recently carried 
out in [34l using a different and slightly smaller parameteriza- 
tion than ours. It is claimed in Ref. 13411 that the addition of 
higher-order spectral parameters does not introduce new param- 
eter degeneracies, and so the results of a Bayesian analysis may 
be expected to yield similar results to those obtained here. 



Physically, e is the equation-of-state parameter of the 
cosmological fluid, and from Eq. ((2]) it follows that the 
condition for inflation, a > 0, requires that e < 1. 

Starting with the equation-of-state parameter, it is 
possible to define an infinite hierarchy of parameters [47j | 
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by taking successive derivatives of the Hubble parameter, 
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In what follows, we will refer to Eqs. {Zj and jS]) col- 
lectively as the Hubble flow parameters. These are of- 
ten referred to as slow-roll parameters in the literature, 
however, they are defined here without any assumption 
of slow-roll. For this analysis it is most convenient to 
use the number of e-folds before the end of inflation, N, 
as our time variable. The scale factor may be written 
a oc e N and from Eq. (fT]), 



dN 



-Hdt=*dS^ "* 



TOPI yfttf) 



(9) 



Making use of this relation, we take successive derivatives 
of the flow parameters with respect to N, generating an 
infinite set of differential equations (3l| , 



dH 

dN 
de 
dN 

(1(7 

dN 
d(?\ H ) 
dN 



eH, 

e(a + 2e), 
-5ecr- 12e 2 



(10) 
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1 



-a+(l-2)e 



where a = 2ij — 4e. In practice, this system is truncated 
at some finite order M by requiring that m+1 Xh = 0. 
This system can then be solved numerically by specify- 
ing the initial conditions of the parameters e, a, . . . , \h 
at some arbitrary time, iVj. Although the system is trun- 
cated at finite order, this results in an exact solution for 
the background evolution of an FRW universe dominated 
by a single scalar field. This is due to the form of the 
flow equations, where it can be seen that the truncation 
ai+1 \h = ensures that all higher-order parameters van- 
ish for all time. In Ref. 31|, the initial conditions were 
drawn randomly from the ranges 



e e [0,0.8] 
a 6 [-0.5,0.5] 

c2 



(11) 



A/+1 A 



G [-0.05,0.05] 
r £ [-0.005,0.005] 

h 6 0, 



although other choices are possible. The system is then 
evolved forward in time until either inflation ends (e > 
1), or the system reaches a late-time asymptote, e — > 0. 
The latter possibility arises in models in which the field 
evolves to a point of nonzero vacuum energy, leading to 
eternal inflation. In such cases, inflation must end via the 
action of an auxiliary field, such as in hybrid models [48| . 
Such models gcncrically predict scalar spectra with n s > 
1 and ne glig ible running, in conflict with the WMAP3 
data-set [22||. We therefore focus solely on models for 
which inflation ends through a failure of slow-roll, e > 1, 
originally termed nontrivial models [3~lj . 

Once a nontrivial model is found, the flow equations 
are integrated backwards in time from N = to N = 
A r G f )s , where N t, s is drawn randomly from the range 
[40, 60] . The solution to the flow equations then comprise 
the full time evolution of the parameters e, cr, . . . , M Xh 
from N = — > N i> s . It is then possible to reconstruct 
the inflaton potential using the Hamilton- Jacobi equa- 
tion, Eq. (U) The values of the flow parameters at 
N i, s can then be used to calculate observables at this 
point via the relations (3r| 



16e[l - C(a + 2e)], 
1 + a - (5 - 3C)e 2 

-i(3-5Cye + i(3-CX 2 , 



(12) 



dn s 
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dn s 
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where C = 4(ln2 + 7) — 5, 7 ~ 0.577, and a denotes the 
running of the scalar spectral index, n s . By adopting 
a Monte Carlo approach, large numbers of models can 
be generated and their observable predictions compared 
with current observational data. This serves as a means 
for constructing classes of models that satisfy certain cri- 
teria, as well as gaining insight into the generic features 
of the inflationary parameter space. Models generated 
stochastically using the above method are observed to 
cluster strongly in the n s — r and a— n s planes [30l. [3ll. |4(| . 
It is difficult, however, to attribute any rigorous statis- 
tical meaning to such clustering, since there is no well- 
defined measure on the parameter space. 

The main drawback of the flow method is that while 
generating exact inflationary evolutions, it must rely on 
approximations when calculating the observable predic- 
tions of these models. These approximations are made 
as two separate series truncations. The first approxima- 
tion has to do with the parameterization of the power 
spectrum. When quoting observables in terms of r, n s 
and a, one is considering a truncated Taylor expansion 
of ln(P(fc)) in ln(fc), 



In 



P(k) 
P(fco) 



(n s - l)ln 
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. (13) 

where k is the comoving wavenumbcr. Clearly, this pa- 
rameterization is only useful if the higher-order terms are 
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small across the range of scales for which it is expected 
to hold. The second approximation results from the use 
of Eq. (fT2|) to connect these spectral parameters to the 
inflationary flow parameters. These expressions are only 
accurate to order C(e 2 ) in the flow parameters. While 
the evolution of the lowest-order parameters appearing 
explicitly in Eq. (|12|) are determined by the full set of 
flow parameters out to order M , the use of these expres- 
sions still requires that these higher-order parameters be 
negligible on observable scales. While, in principle, Eq. 
(fT2")) can be extended to arbitrary order in slow-roll, the 
expressions quickly become very algebraically complex 
[33, HI]. Furthermore, this approach evaluates observ- 
ables at 'horizon crossing', a technique that is not always 
applicable [H, M, ED] ■ 

These approximations arc perfectly valid when con- 
sidering inflation models that satisfy the slow-roll crite- 
ria, e, 77 <C 1. In this case Eq. (fT2")l is highly accurate 
and higher-order spectral parameters are negligible, val- 
idating the use of Eq. (fT5|) . However, if one wishes to 
use flow methods to explore regions of the inflationary 
parameter not well approximated by slow-roll, this ap- 
proach must be adapted. As pointed out by Easther and 
Peiris 42j, if one imposes an inflationary prior on the 
underlying model space, then there is no need to cal- 
culate observables in terms of the spectral parameters. 
The flow parameters themselves fully specify the physics 
of the inflationary epoch and one avoids introducing fur- 
ther, unnecessary errors. We adopt this philosophy in 
what follows, doing away completely with the parame- 
terization Eq. (fT5|) and the specification of observables 
Eq. (HU). 

We develop an improved method of reconstruction by 
combining the flow formalism with a numerical evalua- 
tion of the mode equation of inflationary perturbations. 
By solving the mode equation for multiple fc-values, it 
is possible to reconstruct the primordial power spectra 
without recourse to spectral parameters or the slow-roll 
approximation. This leads to a consistent match be- 
tween the scalar field potential and associated perturba- 
tion spectra, to the same degree of accuracy. In addition, 
it allows for an investigation of regions of the inflationary 
parameter space that lie far from slow-roll, leading to a 
more robust reconstruction program. (For other efforts 
at numerical evaluation of inflationary perturbations, see 
Refs. [I| HQ 53.) 



III. CALCULATING THE PRIMORDIAL 
POWER SPECTRUM 

During inflation, vacuum fluctuations of the inflaton 
field are rcdshiftcd to supcrhorizon scales by the rapidly 
expanding spacetime where they become classical curva- 
ture perturbations. The inflaton field couples at linear 
order to the scalar metric perturbation, 5g s „ v , which may 
be defined in terms of four scalar functions, 



2{H L 5 lJ +d l d J H T )- 



(14) 



If one works in comoving gauge, in which the spatial slices 
are orthogonal to the worldlines of comoving observers, 
the intrinsic curvature perturbation of the spacelike hy- 
persurface is [4{J [5(| 



(15) 



In comoving gauge, one has the additional freedom of 
requiring that the spatial part of the metric perturbation 
be isotropic, Ht = 0, so that 11 = Hl- The coupling 
between inflaton and metric perturbations motivates the 
use of the gauge invariant Mukhanov potential [5l|, [H, 

in, 



aS(j) 



H 



-Mi 



(16) 



where a is the scale factor of the universe, H is the Hubble 
parameter and S(f> is the scalar field fluctuation. Primes 
denote derivatives with respect to conformal time, r. On 
comoving hypersurfaces, S(f> = 0, and the Mukhanov po- 
tential is related to the comoving curvature perturbation, 



Tl = 



(17) 



where z = <p' / H . The two-point correlation function can 
be written in terms of the power spectrum, Pn (k) , 

(ftkftkO = ^-Pn(k)S(k - k'). (18) 

The power spectrum of the curvature perturbation may 
then be written, 



Pn(k) 



2^2 



Uk_ 

z 



(19) 



where the Uk are the Fourier modes of the gauge- invariant 
potential satisfying the equation of motion, 



k 2 



u k = 0. 



(20) 



Tensor metric perturbations, <5.gJ„, are also excited 
during inflation, leading to a large-scale gravitational 
wave background. The metric perturbation is purely spa- 
tial and may be described by a single function, 6gfj = hij . 
The dynamics of these perturbations are determined by 
the linearized Einstein equations, which follow from the 
action 1541. 



c _ m Pl 
bh - 64^ 



(21) 



mi 



2A 



Because gravitational waves are both transverse and 
traceless (i.e. h l i = = 0), they may be decomposed 
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into two independent polarization modes, denoted + and 
x . The Fourier decomposition may then be written 



(22) 



d 3 k 



[/ lk , + (T)e+(k)+/ lkiX (r)eg.(k)] e< kx , 



(2tt) 3 / 2 

where we have used the fact that h\ = /i-k which fol- 
lows from the condition e^k) = e,j(— k). The power 
spectrum of tensor fluctuations can then be be written 

^( fc ) = ^((IK+| 2 ) + (IVx| 2 )). (23) 



Following the field redefinition, 



i 1 

—n 



32tt 



+ ,X ! 



(24) 



the equation of motion that follows from Eq. (|2lj) for 
each polarization mode becomes that of a canonically 
normalized masslcss scalar field in an FRW background, 



v'l + ( k 2 - — ) 



(25) 



where the polarization indices have been suppressed. 
Here v k is the quantum mode associated with positive- 
frequency excitations of the metric perturbation. In 
terms of the new field, the power spectrum Eq. 
becomes 



Ph(k) 



32k 3 v k \ 2 



(26) 



Although written in terms of a single solution of Eq. (f2"5)l , 
this expression includes the contributions from both po- 
larization modes. 

Our goal is to combine Eqs. (|20[) and (|23|) with the flow 
equations, Eq. (| 10|) . Then, for any initial point in the 
flow space, we will be able to determine the full evolution 
of the modes u k and v k together with the background 
evolution. Since the most convenient time variable for 
use with the flow formalism is N, the number of e- folds 
before the end of inflation, it is necessary to recast Eqs. 
(|20|) and |25|) in terms of N. These expressions become 



d 2 u k ^,du k 



and 

d 2 v k ^dv k 



u k = 0, 



(27) 



an ' 



v k =Q, (28) 



where the function F(e, a, £ 2 ) of Eq. |27|) is defined as 
F(e, a,e)-2(l-2e-^-e 2 + '-a 2 + M ■ (29) 



The full system of differential equations is formed by Eqs. 
([27)) and (J25J) together with the flow equations, Eq. (TTU]) . 
The standard choice of initial conditions for the mode 
function is that defined by the Bunch-Davies vacuum, 



u k (-kr -> oo) = \j —< 



— iki 



(30) 



When solving the mode equations numerically we cannot 
use this exact condition for two reasons. First, we cannot 
set this condition in the infinite past, but must impose 
it at a finite time. Imposing the limit Eq. (|5D|) at a fi- 
nite time, i.e. finite length scale, results in modulations 
of the power spectrum akin to those arising from trans- 
planckian effects [5o| . This effect can be minimized by 
initializing the mode functions at sufficiently early times 
(small length scales). Second, we must write Eq. ([30)) 
in terms of t(N) for use with Eqs (J27J) and l[2"5p. and 
this function is not known in general. From the relation 
dN = —aHdr, we obtain 



If--*- 1 '- 



(31) 



where y = k/aH is the ratio of the Hubble radius to the 
proper wavelength of the fluctuation and is a function of 
N. If e(y) is approximately constant, this equation can 
be integrated to give r(y). From the equation 



de 
dy 



1 



de 



2/(1 -e) dN 1 



(32) 



we see that e(y) = const, if y is taken sufficiently large. 
By taking y large, we are also ensuring that the modes are 
initialized in the short- wavelength limit. This motivates 
the use of the approximate initial conditions, 



du k 
dN 



2k 

T 



-iVi/(l-u) 



(33) 



For a choice of initial flow parameters e^, <7j, Xh%, 
Eq. (fTl"]) . we set the initial condition for each fc-mode at 
— e,) = 100. This proves to be sufficiently large to 
ensure the accuracy of the conditions Eq. . 

Since there are two complex solutions to each of Eqs. 
(|27p and l|28p. rather than work with complex coefficients, 
we define the orthogonal solution basis, 



Uk,l 
U k ,2 



u k + u* k 
2 

u k - u* k 
2i ' 



(34) 



Each mode is evolved from N(yi/(l—€i) = 100) to N = 0, 
the end of inflation. The amplitude of the power spec- 
trum for each fc-mode is then obtained by evaluating Eqs. 
(TiU)) and (|26[) at this time, when all modes have attained 
the long- wavelength limit. 
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IV. SPECTRUM RESOLUTION WITH WMAP3 

In this section we investigate the ability of the latest 
WMAP data-set to accurately resolve the form of the 
power spectrum. We approach this problem by consider- 
ing a 7-parameter hducial best-fit model parameterized 
by inflationary and non-inflationary degrees of freedom. 
We suppose that this model accurately describes the uni- 
verse and search for alternative models which arc statis- 
tically indistinguishable from it. The non- inflationary 
parameters are taken to be the baryon and CDM densi- 
ties, fifc/i 2 and fl c h 2 , the Hubble parameter, h, the op- 
tical depth to rcionization, t, and the overall spectrum 
normalization, A(k = 0.002 /iMpc -1 ). 2 The inflationary 
parameters are the scalar spectral index, n s , and the ten- 
sor contribution, r. We consider purely adiabatic initial 
perturbations and assume spatial flatness. The tensor 
spectral index, rit, is assumed to satisfy the inflationary 
consistency condition, n t = —r/8, and docs not represent 
an additional free parameter. 

Since we only vary the form of the spectrum in this 
analysis, we fix the non-inflationary parameters at their 
best-fit values. We then replace the power-law parame- 
terization (n s , r) with tensor and scalar power spectra 
generated with the Monte Carlo method. We work to 
6 t,l -order in the flow space, i.e. 1 \h = for i > 6. We 
use a version of CAMB [35[ modified to accept arbitrary 
power spectra as input to generate the associated Cg- 
spectra. It is then possible to calculate the model's effec- 
tive chi-square, Xeff = — 21n£, using the WMAP3 likeli- 
hood software available at the LAMBDA website 3 . One can 
then collect models that satisfy specific likelihood criteria 
relative to the best-fit model. 

How does one compare the statistical significance of 
two different models? If the two models have the same 
parameterization, then the confidence limits obtained 
from a maximum likelihood analysis of the parameter 
space suffice. However, in this study, we wish to compare 
a model comprising 2 spectral parameters with a model 
comprising 6 flow parameters. A simple statistic used 
in model selection analyses is the likelihood ratio test, 
21n£ s i m pi e /£ comp i ex , where the simple model is so named 
because it contains fewer free parameters than the more 
complex model. This statistic is approximately xt dis- 
tributed with v degrees of degrees of freedom equal to 
the difference in the number of free parameters between 
the complex and the simple model. The significance level 
associated with the value of the xt can then be used to 
quantify the relative goodness-of-fit between these mod- 
els. Use of this test requires that the two models be 
nested - that the complex model is formed by adding 



2 Since the overall spectrum amplitude is not fixed by the in- 
flationary model in the flow formalism, it is considered a non- 
inflationary parameter in this analysis. 

3 http://lambda.gsfc.nasa.gov/ 



parameters to the base model. The space of flow pa- 
rameters and spectral parameters are at best only ap- 
proximately nested: there is no way to map the flow 
parameters (e, cr, £ 2 , 3 A#, 4 Xh, 5 A#) to a finite number of 
spectral parameters. This is because there are infinitely 
many higher-order spectral parameters that are functions 
of these 6 flow parameters that in general will not vanish 
under the mapping. We say approximate because if these 
parameters are very small, they might be neglected. A 
further downside is that it is not clear how to tell from 
this statistic if two models are of comparable significance. 

Another method often used in model selection analyses 
is the Bayesian information criterion (BIC) [57j |. 

BIC = -21n£ + fctoJV, (35) 

where C is the maximum likelihood of the model, k the 
number of free parameters and N the number of data 
points. The BIC is an approximation of the Bayesian 
evidence, which is the integral of the likelihood function 
over the full parameter space. The BIC penalizes overpa- 
rameterized models that don't provide significantly bet- 
ter fits to the data. This approach is also ill-suited for 
our purposes, because it assumes that all k parameters 
are well measured by the data. In fact, it is the oppo- 
site of this case that motivates this study: current data 
is not good enough to resolve higher order terms in the 
spectral decomposition, and we seek to determine what 
these resolution limits are. 

In what follows, we instead consider the p- value calcu- 
lated from a model's chi-square per degrees of freedom as 
a measure of goodness-of-fit. This method is applicable 
so long as the likelihoods are approximately Gaussian. 
The p-value of a proposed model with a given \ 2 is 

/>OC 

p v = / P v {y)dy, (36) 

where P v is the chi-square probability distribution func- 
tion with v degrees of freedom (d.o.f.), 

Pu{y) = rf^ 72 " 1 ^ 2 - (37) 

For a given significance level, a, if p < a then the pro- 
posed model may be rejected at the 1 — a confidence level. 
More precisely, the p-value is the probability that we ob- 
tain a particular \ 2 due to chance alone. It can there- 
fore be interpreted as the probability of falsely rejecting 
a correct model of the universe. To determine relative 
goodness-of-fit between two models, one simply compares 
their p- values. In order to be as conservative as possible, 
we consider trial models lying within |Ap| < 0.01 of the 
best- fit model to be effectively indistinguishable from it. 
The spectra generated by G^-order Monte Carlo are de- 
scribed by 6 free parameters, corresponding to the flow 
parameters e, a, 5 \h- Therefore, the trial spectra 
contain 4 more parameters than the best-fit power-law 
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FIG. 1: A representative sampling of power spectra lying 
within |Ap| < 0.01 of the best-fit model. The spectrum of 
the best-fit model is represented by the heavy dashed line. 
The spectra are colored-coded according to their predicted 
values for r at k = 0.002 ftMpc -1 : red yield r ~ O(10 _1 ) and 
black yield r < O(10 _1 ). 



model. The difference in p-value between a trial model 
with v — 4 d.o.f. and the best-fit model with v d.o.f., 



|Ap| = <0.01, 



(38) 



corresponds to the likelihood spread 
- 8 x 10- 4 < ^,J{y - 4) - Xbest-fit/^ £ 10" 3 - (39) 

The chi-square per degrees of freedom of the best-fit 
model is xL t _ fi > = 1-0210. 

We begin with a best-fit model with n s = 0.969, 
A = 2.30 x 10 -9 and a tensor-to-scalar ratio, r = 0.0346. 
In Figure [1] we present a sample of power spectra ly- 
ing within \Ap\ < 0.01 of this model, and in Figure [2] 
some reconstructed potentials. In order to identify spec- 
tra exhibiting the strongest deviation from power-law, 
we initially imposed a non-slow-roll prior on the model 
space. In order for a model to be considered for full power 
spectrum integration, we require that the spectral index 
calculated to 2™ d -order in slow-roll Eq. (fT2|) differ from 
the 3 rd -order result [H| by at least 1%. This search re- 
sulted in the grouping of red-colored spectra exhibiting 
a large suppression of power on large scales in Figure Q] 
We later relaxed this prior to obtain the black-colored 
spectra, allowing us to form a degeneracy envelope. 

Remarkably, the red spectra remain equally good fits 
to the data as the fiducial power-law even when SDSS 
data [HI] is included. For the case of power-law spec- 
tra, the latest analyses report r .oo2 < 0.30 at 95%- 
confidence when SDSS is combined with WMAP3, while 
T-o.002 < 0.65 at 95%-confidence with WMAP3 alone [12]. 
The error bars open up considerably when running is al- 
lowed: ro.002 < 0.38 at 95%-confidence SDSS+WMAP3, 



FIG. 2: A representative sampling of reconstructed potentials 
lying within |Ap| < 0.01 of the best-fit model. The potentials 
are colored-coded according to their predicted values for r at 
k = 0.002 /iMpc" 1 : red yield r ~ O(10 _1 ) and black yield 
r < O(10 _1 ). The potentials have been given a common, 
arbitrary normalization at = 0, when scales corresponding 
to the quadrupole exit the horizon. 



versus ro.002 < 1-1 with WMAP3 alone. Spectra with 
larger ro.002 can be accommodated by the data if there is 
suppressed scalar power on these scales. As r decreases 
on smaller scales, the scalar spectrum must run to larger 
amplitudes in order to maintain the correct amount of 
overall power. The SDSS data provides an accurate 
measurement of the matter power spectrum on scales 
0.01 MVIpc" 1 < k < 0.3 Mvlpc" 1 , and while not prob- 
ing scales k ~ 0.002 /iMpc -1 directly, if there is a strong 
negative constant running of the spectrum on these scales 
it may lead to a measurable change in power on the inter- 
mediate scales that are directly probed by SDSS. In par- 
ticular, if the running is too large at k ~ 0.002 /iMpc -1 , 
there will be a loss of power on smaller scales as the spec- 
trum dips blue again. This results in the much tight- 
ened bound on r. While the red spectra in Figure [T] 
arc strongly running on large scales, permitting a size- 
able tensor-to-scalar ratio, this running is not constant 
and is unsubstantial on small scales. Indeed, when the 
red spectra are plotted in the n s -r plane, we find that 
they all lie outside the 1-a contour of the SDSS+WMAP 
marginalized likelihood, Figure [3] Therefore, by consid- 
ering more general spectra, the error bars inevitably get 
much larger. The spectral index of the non-power-law 
spectra is calculated from the slope of the spectrum at 
k = 0.002 /iMpc" 1 , 



dlnP(fc) 



d\nk 



(40) 



fc=0.002 



This is a measure of the local value of n s around k = 
0.002 /iMpc" 1 and it should be emphasized that the large 



8 



0.5 - 




n 

s 



FIG. 3: Constraints on the n B -r plane from combined SDSS 
and WMAP data for spectra with constant running. The 
inner and outer contours mark the 68%- and 95%-confidence 
intervals, respectively. The red points denote the values of n a 
and r of the red spectra in Figure [1] The blue square is that 
of the fiducial power-law model. 




FIG. 4: The best-fit inflaton potential (blue dashed line) as 
compared to a model lying within Ap < 0.01 of this model. 
The potential energy of the alternative model is given by the 
red solid line, and the kinetic energy is given by the red dot- 
ted line. The vertical line denotes the time when scales cor- 
responding to the quadrupole leave the horizon. 



uncertainty is due to the generality of the spectral shapes. 
However, the uncertainty in the local value of n s is dras- 
tically reduced at smaller scales, as indicated in Figure 
[TJ Recently, Cortes et dl. [58| determined that the effect 
of the degeneracy in the n s -r plane can be minimized by 
quoting observables at k = 0.017 /iMpc -1 . 



The statistical degeneracy of the power spectra exhib- 
ited in Figure [1] translates directly to a degeneracy in the 
space of inflation models, Figuref^] The red-colored infla- 
tion models that give rise to the red-colored non-power- 
law spectra are characterized by an initially fast-rolling 
inflaton that has failed to reach the slow-roll attractor 
by the time observable scales exit the horizon, similar 
to models proposed by Contaldi, et al. 59]. In Figure 
[4] we plot one such potential. The kinetic energy is ini- 
tially large and monotonically decreasing with respect to 
time. From the equation of motion Eq. ([3]), it is evi- 
dent that if the field velocity is initially large, the Hub- 
ble drag term 3H(f> dominates the evolution, slowing the 
field. The field decelerates until the drag term becomes 
subdominant near the end of inflation. In contrast, the 
kinetic energy of the best-fit power-law model (not in- 
cluded in the figure) is 0(1O~ 15 ) and is a monotonically 
increasing function of time. 

We conclude that current CMB data does not allow for 
a precise reconstruction of the inflationary power spec- 
trum and the corresponding inflaton dynamics. In partic- 
ular, we find a statistical degeneracy amongst slow-roll 
and non-slow-roll models. The span of the degeneracy 
envelope in Figure [1] is largely shaped by our inability to 
constrain r, with models exhibiting the strongest devia- 
tion from power-law also predicting values of r an order of 
magnitude larger than the fiducial model. Future CMB 
missions 0, [2J] with the ability to detect the B-mode 
polarization signal characteristic of gravitational radia- 
tion will put tighter constraints on r [60l . [6lT | , reducing or 
breaking the r-n s degeneracy and significantly improving 
primordial power spectrum resolution. 



V. CONCLUSIONS 

We have investigated the ability of current data to con- 
strain the form of the primordial power spectrum. The 
ability of the data to resolve the power spectrum can be 
tested by how well it singles out a best-fit model. Start- 
ing with the latest WMAP release, we generate a best-fit 
power law model consistent with the simplest models of 
slow-roll inflation. We test a wide variety of more com- 
plicated power spectra, and identify those that provide 
equally good fits to the data relative to the best-fit model. 
If the alternative spectra arise in an inflationary context, 
then we can hope to identify both slow-roll and non-slow- 
roll models that are equally consistent with current data. 

In order to test a wide array of inflationary power spec- 
tra, we turn to Monte Carlo reconstruction. We com- 
bine the flow formalism, which is a method of stochastic 
model generation, with a numerical integration of the 
mode equations of quantum fluctuations. This allows us 
to handle inflation models that yield spectra that are 
not well described by the standard spectral parameter- 
ization. For each spectrum thus generated, we fix the 
non-spectral parameters at their best-fit values and cal- 
culate the likelihood of the model. By fixing the non- 
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spectral parameters, we are only sampling a subset of 
possible spectra that might be degenerate with the best- 
fit power-law model, making this approach conservative. 
We determine the statistical significance of each model 
by obtaining the p-value calculated from the Xg^j/d.o.f., 
Eq. (J36|). 

We generate an ensemble of power spectra to 6 th - 
order in the flow space and select only those lying within 
Ap < 0.01 of the best-fit model. A sampling of power 
spectra meeting this criterion are shown in Figure [T] The 
current CMB data provided by WMAP3 [H| only reliably 
constrains the form of the power spectrum on intermedi- 
ate scales, 0.01 /iMpc -1 < k < 0.1 /iMpc -1 , with much 
variation on larger spatial scales where cosmic-variance 
is the dominant source of error. By doing away with the 
spectral parameterization Eq p3[) . we also free ourselves 
from the constraints imposed on the tensor-to-scalar ra- 
tio by SDSS data, the tightest bounds affecting models 
with a constant running of the spectral index. The flow 
method allows us to generate models with varied forms, 
including models with non-constant running. Such mod- 
els can support a relatively large value of ro.002 by hav- 
ing a large running on these scales that turns off on the 
intermediate scales probed by SDSS. For example, we 
find spectra with ro.002 ~ 0.5 that yield equally good fits 
to the WMAP+SDSS data as the best-fit power-law, a 
value which is well outside the r < 0.38 95%-confidcncc 
limit for spectra with running. By considering more gen- 
eral power spectra, one therefore opens up regions of pa- 
rameter space excluded in simpler models. The inflation 
models responsible for generating the strongly running 
spectra in this study are characterized by an initially 
fast-rolling inflaton. Observable scales exit the horizon 
before the field slows to the slow-roll attractor, yielding 
a large amplitude of gravitational waves and significant 



running. We arc able to conclude that while slow-roll 
inflation models yield perhaps the simplest explanation 
for the origin of large scale structure, fast-rolling infla- 
ton fields are equally suitable candidates from a strictly 
data-driven standpoint. Another interesting approach is 
to expand the inflationary model space beyond the sim- 
plest class of canonical single field models. An example 
is the DBI inflation scenario [62|, [HI analyzed using flow 
techniques by Peiris et al. [64[ , who find that such models 
are strongly constrained by the existing data. 

Finally, we draw attention to the work of Lesgourgues 
and Valkenburg, Ref. [34|. In that analysis, the authors 
consider an enlarged inflationary parameter space con- 
sisting of the coefficients of a Taylor expanded inflaton 
potential. They perform a Bayesian parameter estima- 
tion analysis and obtain confidence intervals on an en- 
semble of power spectra and their associated inflaton po- 
tentials. This choice of parameterization excludes models 
that deviate strongly from slow-roll, in particular, mod- 
els of the type found in this analysis. Aside from this 
difference, we find good agreement between our results, 
which arc based on inferential statistics that are frequen- 
tist in nature, and the results of Ref. [HI, derived using 
strictly Bayesian methods. The bottom line is that cur- 
rent CMB data does not reliably constrain the form of 
the power spectrum, and this conclusion can be reached 
from either a Bayesian or a frequcntist approach. 
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